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ABSTRACT 

By means of a Monte Carlo code that accounts for Compton scattering and photoabsorption followed by recom- 
bination, we have investigated the radiation transfer of Lya, Hea, and recombination continua photons of H- and 
He-like C, N, O, and Ne produced in the photoionized atmosphere of a relativistic black hole accretion disk. We 
find that photoelectric opacity causes significant attenuation of photons with energies above the O VIII K-edge; 
that the conversion efficiencies of these photons into lower-energy lines and recombination continua are high; and 
that accounting for this reprocessing significantly (by factors of 21% to 105%) increases the flux of the Lya and 
Hea emission lines of H- and He-like C and O escaping the disk atmosphere. 

Subject headings: accretion, accretion disks — black hole physics — galaxies: individual (MCG —6-30-15) — 
galaxies: Seyfert — radiative transfer — X-rays: galaxies 



1. INTRODUCTION 

A controversy currently exists concerning the interpretation 
of the broad spectral features observed in the soft X-ray spectra 
of type 1 Seyfert galaxies. Lee et al. (2001) and Turner et al. 
(2003) interprete the features seen in the Chandra High Energy 
Transmission Grating (HETG) and XMM-Newton Reflection 
Grating Spectrometer (RGS) spectra of MCG -6-30-15 as ab- 
sorption edges produced by a dusty, partially ionized ("warm") 
absorber, while Branduardi-Raymont et al. (2001), Mason et al. 
(2003), and Sako et al. (2003) interpret the features seen in the 
XMM-Newton RGS spectra of MCG -6-30-15 and Mrk 766 
as relativistically broadened and gravitationally redshifted Lya 
emission lines of H-like C, N, and O arising in the accretion 
disk surrounding a maximally spinning Kerr black hole. While 
warm absorbers are a common feature of Seyfert 1 galaxies, the 
relativistic line interpretation is appealing because it produces 
model parameters (disk inclination, emissivity index, and inner 
radius) that are not inconsistent with those derived for the Fe 
Ka line. 

To be consistent with observations, proponents of the rela- 
tivistic line interpretation must explain the apparent absence of 
the Ly/? lines and radiative recombination continua (RRCs) of 
C, N, and O, as well as emission by higher-Z ions. Branduardi- 
Raymont et al. (2001) and Sako et al. (2003) propose that this 
is accomplished by radiation transfer effects, specifically that 
the higher-order Lyman lines are suppressed by line opacity, 
while the C, N, and O RRCs and the higher-Z emission fea- 
tures are suppressed by photoelectric opacity. The production 
of C, N, and O emission lines in photoionized constant-density 
and hydrostatic disk models has been studied quantitatively by 
Ballantyne, Ross, & Fabian (2002, hereafter BRF). These cal- 
culations do not, however, account for the reprocessing of soft 
X-ray photons into lower-energy lines and RRCs. In a mod- 
erately ionized solar abundance plasma, the continuum opacity 
is dominated by C from 0.39 to 0.74 keV and O from 0.74 to 



2.4 keV (Fig. 1). If, for example, the medium is optically thick 
to photoelectric absorption by 0^+, a fraction of the photons 
with energies greater than that of the O VIII K-edge will be ab- 
sorbed, and a fraction of those will be reradiated as O VIII Lya 
photons, thereby increasing the equivalent width of the line. To 
explore such processes quantitatively, we have combined a disk 
atmosphere code with a Monte Carlo radiation transfer code to 
calculate the emergent spectrum of the irradiated atmosphere of 
a relativistic black hole accretion disk. In this first communica- 
tion, we concentrate on the physics of the radiation transfer of 
Lya, Hea, and recombination continua photons of H- and He- 
like C, N, O, and Ne, leaving to a future work the calculation of 
equivalent widths and detailed comparisons with data. 

2. DISK ATMOSPHERE AND RADIATION TRANSFER CODES 

The disk atmosphere code calculates the structure of the at- 
mosphere of a black hole accretion disk photoionized by an ex- 
ternal source of nonthermal X-rays. The code is adapted from 
that used by Jimenez-Garate, Raymond, & Liedahl (2002) to 
model the accretion disk atmospheres of X-ray binaries. It as- 
sumes axial symmetry and solves the hydrostatic balance, pho- 
toionization equilibrium, and radiation transfer equations on an 
adaptive mesh in the slab approximation in a succession of de- 
coupled disk annuli using the Raymond (1993) photoionization 
plasma code, which includes the X-ray emissivities and opaci- 
ties of all ions of the 12 most abundant elements. The disk is as- 
sumed to have solar abundances, to be gas-pressure dominated, 
and to have a structure given by Shakura & Sunyaev (1973) 
with the scale radius r» = = GAI / (?, modified by the Page 
& Thorne (1974) and Riffert & Herold (1995) general relativis- 
tic corrections to the vertical structure. We assume that the in- 
cident angle of the nonthermal radiation i = 55°, to simulate 
isotropic illumination, and that the ratio of the incident nonther- 
mal flux to the disk thermal flux / = / Fxdv j J F^iskdv = 1. 
The lower boundary conditions of the calculation are the disk 
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temperature [TdiskC^") oc r^^/'^] and pressure, while the up- 
per boundary condition is the Compton temperature, which is 
reduced to Tc ~ 10 K by the thermal disk emission that irra- 
diates the atmosphere from below. 

As with LMXB disks, the structure of the atmosphere of an 
AGN disk is sensitive to the assumptions about the behavior of 
photoionized plasma on the unstable branches of the thermal 
equilibrium S-shaped curve (Krolik, McKee, & Tarter 1981). 
To maximize the production of soft X-ray lines from H-like 
ions, we assume thermal equilibrium solutions on the upper 
branch of the S-shaped curve. Furthermore, we ignore the ef- 
fects of the radiation pressure in the hydrostatic balance equa- 
tion, since it reduces the thickness of the atmosphere and the 
production of soft X-ray lines, and since in a real disk dynam- 
ical and magnetic pressures may counteract such a reduction 
(Miller & Stone 2000). We assume that the nonthermal con- 
tinuum has a power-law shape with a photon index F = 2.1 
with a high-energy exponential cutoff ii'high = 150 keV and a 
low-energy exponential cutoff i?iow = ^^disk; the black hole 
has mass M — 10^ M© and spin a = 0.998, so the radius 
of the innermost stable circular orbit nsco = 1-23 rg (Thorne 



1974); the accretion rate M = 10^'^ g = 0.2 MEdd; and the 
viscosity parameter a = 0.01. The model has 499 z x 35 r = 
17.5 X 10"^ cells, with the annular radii r = 1600 rg, with 
k = 2,3,4,. . . ,36. For additional details about the disk atmo- 
sphere model, see Jimenez-Garate et al. (2004). 

While the radiation transfer equation used in the disk atmo- 
sphere code is sufficiently accurate to calculate the disk atmo- 
sphere structure, it cannot be used to calculate its spectrum. To 
solve this problem, we are developing a Monte Carlo code that 
accounts for radiation transfer in the atmosphere of a relativistic 
black hole accretion disk. On a microphysical scale, we assume 
that photons are subject to Compton and photoelectric opacity. 
For a photon of energy E scattering off an electron drawn from 
a population with an isotropic Maxwellian velocity distribution 
with temperature kT, the Compton cross section 

r E \ 26/ S \2/, kT 
ac=OT 1-2 J +^ ^ 1 + 5 

where ut is the Thomson cross section, which is valid for E <C 
meC^ and kT < m^c? (Rybicki & Lightman 1979). For a pho- 
ton with energy E and direction vector p scattering off an elec- 
tron with velocity 3c, 7 = (1 — /3^)~^/^, and direction vector e, 
the scattered photon energy E' = 11' (0) /c and direction vector 
p' = n'(l:3)/n'(0). In these expressions, the scattered photon 
4-momentum H' = R''^ l^-^R-^S{l\e)R^kRJ\, where i?e 
and R^ are the electron and photon rotation matrices, respec- 
tively; A is the Lorentz transformation matrix; and S{A6) = 
QR{A6) is the scattering matrix. In this expression, Q = 
1/ [1 +^E ( 1 - /3 e -p) (1 - cos A6) / rueC^] and the scattering an- 
gle AO is given by the Klein-Nishina formula for the differential 
Compton scattering cross section. The photoelectric opacity is 
calculated at each energy for 446 subshells of 140 ions of the 
12 most abundant elements of a solar abundance plasma using 
the partial photoionization cross sections of Vemer & Yakovlev 
(1995). Following photoabsorption by K-shell ions, we gen- 
erate RRC and recombination line cascades in a probabiUstic 
manner using the recombination cascade calculations described 
by Sako et al. (1999). The shapes of the RRCs are determined 
by the functional form of the photoionization cross sections and 
the local electron temperature. Recombination emission from 
L-sheU ions is ignored, since it is suppressed by resonant Auger 
destruction (Ross, Fabian, & Brandt 1996). 



On a macrophysical scale, we assume that the velocity of 
the gas is restricted to the azimuthal direction and has a value 
V4, = vk/[^ + {vk/c)^], where vk = y/GMjr is the Keple- 
rian value, as appropriate for a maximally spinning Kerr black 
hole (Bardeen, Press, & Teukolsky 1972); that within each cell 
the gas moves with a uniform velocity; and that photons travel 
in straight Unes. By imposing a limit on the velocity differ- 
ence from one cell to the next, we determine the radial and 
azimuthal dimensions of the cells used in the Monte Carlo cal- 
culation. Specifically, we require that Av < 2 x 10^ cm s^^ 
[set by At; = AAc/A with A = 1.78 A (Fe XXVI Lya) and 
AA = 0.012 A (the Chandra High Energy Grating FWHM)]. 
Extending the grid from r = 1.23 to 600 Vg, the model has 

499 5 X 108 f X 52-1005 ^ = 21.5 x 10'^ cells. 

Photons are propagated in this grid until they are destroyed or 
escape. From the "event list" of escaping photons, we can cal- 
culate the emergent spectrum in the fluid frame, in the local in- 
ertial frame (accounting for Lorentz boosts), and in the inertial 
frame at infinity (accounting for gravitational redshifts). How- 
ever, because we do not yet follow photons along geodesies, 
we are not yet able to directly calculate the observed spectrum 
at infinity at a given disk inchnation. To solve this problem, 
we extracted from the Laor (1991) model in XSPEC the line 
profiles Lk iy) and relative weights for a (5-function spec- 
trum for each disk annulus for an inclination angle i = 35°, 
as typical of Seyfert 1 galaxies. The observed spectrum at in- 
finity at this inclination in the Kerr metric is then S'ivi) = 
J2k'^kJ2j Sk{i'j)Lk{vi,yj), where Sk{i>) are the spectra in 
the fluid frame for each disk annulus. 
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Fig. 1 . — Energy dependence of the optical depth to photoelectric absorption 
of a solar abundance plasma with unit ionization fractions of H- and He-like C, 
N, O, Ne, Mg, and Si. Energies of the Lya (Hea) lines of H-like (He-like) C, 
N, O, Ne, Mg, Si, and S are indicated. Note that C (O) dominates the opacity 
from 0.39 to 0.74 (0.74 to 2.4) keV. 



3. MONTE CARLO SPECTRA 

The Monte Carlo code can be configured to run in a number 
of ways, but for the present investigation we used it to study 
the radiation transfer of recombination photons produced in the 
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disk atmosphere (tt < 10). For each transition, the emissiv- 
ities jz,i = neUpAz fz,i+iC(T^^ photons cm^"^ s^^, where 
rie is the electron density, Hp is the proton density, is the 
abundance of element Z relative to H, fzd+i is the ionization 
fraction of the i + 1 charge state of element Z, and aT''^, 
is the temperature-dependent "per-ion" recombination rate co- 
efficient based on power-law fits to the recombination rates 
over the temperature range of interest. Because the recombi- 
nation luminosity of our disk atmosphere model is dominated 
by low-Z ions, we performed Monte Carlo calculations for the 
n = 2 — » 1 (Lya doublets and Hea triplets) and RRC tran- 
sitions of H- and He-Uke C, N, O, and Ne. In addition, we 
performed calculations for the Lya doublets of H-like Mg, Si, 
S, Ar, Ca, and Fe and for the brightest line of the brightest Fe 
L-shell ion (Fe XVII A17.05 A), but found for the adopted disk 
model that they have only a small effect on the net spectrum. 

Figure 2 shows the emergent spectra resulting from the pro- 
duction in the disk of Lya photons of H-like O and Ne. In the 
upper panel, 80% of the escaping photons are in the O VIII Lya 
line and 9% are in the C VI Lya line. In the lower panel, 42% 
of the escaping photons are in the Ne X Lya Une, 26% are in 
the O vm Lya line, and 5% are in the C VI Lya line. The 
conversion efficiency of Ne X Lya into O VIII Lya is so high 
because the former transition hes just above the O VIII K-edge 
(Fig. 1). In both cases, very Uttle flux escapes in the N recombi- 
nation lines because the C and O opacities dominate over the N 
opacity, and N recombination lines are subsequently converted 
into C recombination Unes. 
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Fig. 2. — Emergent spectra resulting from the production in the dislc of Lya 
photons of H-like O (a) and Ne (i>). Black curves are the spectra in the fluid 
frame, solid gray curves are the spectra at infinity at i = 35°, and the dotted 
gray curves are the spectra ignoring recombination following photoionization. 
Energies of various Lya and Hea lines are indicated. 

Figure 3 shows the emergent spectra resulting from the pro- 
duction in the disk of RRC photons of H-like O and Ne. In 
the upper panel, 39% of the escaping photons are in the O VIII 
RRC, 30% are in the O Vni Lya fine, and 6% are in the C VI 
Lya line. In the lower panel, 52% of the escaping photons are 
in the Ne X RRC, 4% are in the Ne X Lya line, 16% are in 
the O VIII Lya line, and 3% are in the C VI Lya line. The 
O VIII RRC conversion efficiency is so high because the RRC 
photons are emitted just above the O VIII K-edge, where the 
opacity is the highest. The conversion efficiency of Ne X RRC 



photons into the Ne X Lya line is so low because the O opacity 
dominates over the Ne opacity at the Ne X K-edge (Fig. 1), and 
because Ne X Lya photons are subsequently converted into the 
O Vin Lya line with high efficiency (Fig. 2V). Once again, the 
efficiency of conversion into escaping N recombination Unes is 
very low. 
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Fig. 3. — Similar to Figure 2, for the RRC photons of H-like O (a) and Ne 
(i>). Energies of various Lya and Hea lines and the RRCs are indicated. 

As these examples demonstrate, recombination following 
photoionization reprocesses soft X-ray photons into lower- 
energy lines and RRCs, and in the process increases the 
strengths and hence the equivalent widths of the C and O emis- 
sion features. As shown in Table 1, for the adopted disk model, 
the efficiencies of conversion of the Lya, Hea, and RRC pho- 
tons of H- and He-like C, N, O, and Ne into lower-energy lines 
and RRCs range from 20% (for O VIII Lya) to 86% (for Ne IX 
Hea). 

TABLE 1 



Conversion Efficiencies of Escaping Photons 



Transition 


C 


N 


O 


Ne 


H-UkeLya 




36% 


20% 


58% 


He-like Hea 




30% 


39% 


86% 


H-like RRC 


41% 


54% 


63% 


72% 


He-like RRC ... 


. 34% 


33% 


61% 


48% 



To show the net effect of these processes on the observed 
spectrum, we show in Figure 4 the full disk atmosphere spec- 
trum accounting for the n = 2 ^ 1 and RRC transitions of 
H- and He-like C, N, O, and Ne, with relative normaUzations 
determined by the integrated disk atmosphere emissivities. In 
this model, the emission spectrum is dominated by C, N, and 
O, so the inclusion of Ne (and higher-Z ions) makes only a mi- 
nor contribution to the net soft X-ray spectrum. The emergent 
spectrum is shown two ways in Figure 4. In the upper panel, 
recombination following photoionization is ignored, while in 
the lower panel this process is fully accounted for. These spec- 
tra differ in two important ways. First, the lower spectrum is 
brighter than the upper spectrum by 50%: by accounting for re- 
combination following photoionization, significantly more line 
flux escapes the disk atmosphere. Specifically, for the adopted 
disk model, the flux and hence the equivalent widths of the C V 
Hea, C VI Lya, O VII Hea, and O VIII Lya lines increase 
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by 21%, 105%, 28%, and 44%, respectively. Second, the rela- 
tive line strengths change: the C VI Lya to C V Hea line ratio 
increases by 69%, the O VIII Lya to O VII Hea Une ratio in- 
creases by 12%, and the summed C/O Une ratio increases by 
16%. In both cases, the N emission features are weak because 
the N emissivities are relatively low and little O and Ne flux is 
reprocessed into N recombination features (note that the emis- 
sion feature at 0.5 keV in Fig. 4 is a combination of N VII Lya 
and the C VI RRC). 
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Fig. 4. — Emergent .spectra resulting from the production in the disk of Lya, 
Hea, and RRC photons of H- and He-like C, N, 0, and Ne. Panel (a) ignores 
recombination emission following photoionization, while panel (h) accounts 
for this emission. Black curves are the spectra in the fluid frame and gray 
ciu-ves are the spectra at infinity at i = 35°. Energies of various Lya and Hea 
lines and the O RRCs are indicated. 

4. SUMMARY AND CONCLUSIONS 

By means of a Monte Carlo code that accounts for Compton 
scattering and photoabsorption followed by recombination, we 
have investigated the radiation transfer of Lya, Hea, and RRC 
photons of H- and He-like C, N, O, and Ne produced in the 
photoionized atmosphere of a relativistic black hole accretion 
disk. We find that the conversion efficiencies of these photons 
into lower-energy fines and RRCs are higfi, and that accounting 
for this reprocessing significantly (by factors of 21% to 105%) 
increases the flux of the Lya and Hea emission lines of H- and 
He-like C and O escaping the disk atmosphere. These fines 
dominate the observed spectrum (Fig. 4b) because the relative 
emissivities of these ions are high; the relative photoelectric 
opacities of these ions are high (Fig. 1), so photons with en- 
ergies of 0.39 to 2.4 keV are preferentially absorbed by C and 
O and converted into fine and RRC photons (Fig. 2); and RRC 



photons are converted with high efficiencies into lines (Fig. 3). 
The relative strength of the C VI Lya line is particularly high 
in our model because it is fed strongly by aU higher-Z ions. 

Our model spectrum appears to differ from the Chandra 
HETG <mdXMM-Newton RGS spectra of MCG -6-30-15 and 
Mrk 766 in having too much flux in the C V and O VII Hea 
fines and too little flux in the N VII Lya line. The relative 
strength of the C V and O VII Hea lines may not be of concern, 
since it should be possible to weaken these lines by increas- 
ing the mean ionization level of the disk atmosphere: in their 
photoionized constant-density disk models, BRF find that the 
Lya fines dominate the soft X-ray spectrum for ionization pa- 
rameters ^ > 500. In contrast, the relative weakness of the 
N VII Lya line may be of concern. BRF find that the strength 
of this line can be increased simply by increasing the N abun- 
dance. However, we note that if the opacity from 0.55 to 0.67 
keV comes to be dominated by N^+, it and not C will cap- 
ture the O VIII Lya and O VII Hea photons, leading to an in- 
crease in the N VI Hea line strength. Sako (2003) has proposed 
an alternate way to increase the N VII Lya line strength by a 
Bowen-like fluorescence mechanism, wherein O VIII Lya pho- 
tons pump the N VII Lyi^ transition, and the resulting N Vn Ly<^ 
photons are converted into N VII Lya photons via line opacity. 
We are not yet able to test this mechanism quantitatively, but we 
expect that it will not be particularly efficient in an irradiated 
disk atmospfiere, since the O VIII Lya photons are produced in 
a layer where the 0^+ ionization fraction is low and the 
ionization fraction is nearly zero. Only those O VIII Lya pho- 
tons that are directed down into the disk and survive passage 
through the 0^^+ layer will have a chance of being converted 
into N VII Lya photons. 

In its current form, our Monte Carlo code treats in consider- 
able detail many of the physical processes that affect the trans- 
fer of X-ray photons in the photoionized atmosphere of a rel- 
ativistic black hole accretion disk, but it does not yet follow 
photons along geodesies and it does not yet account for line 
opacity. When these processes are included in the code, we 
will be able to directly calculate the disk spectrum observed at 
a given inclination at infinity and study tfie Bowen fluorescence 
mechanism, the radiation transfer of higher-order Lyman lines, 
and many other physical processes that affect the X-ray spectra 
of AGN disks. 
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